program fortran

  implicit none
  integer i,j,k,n
  real*8 size, dx
  real*8 lambda, mu, slip, opening, width, length, strike, dipd, rake, origin(3)
  real*8 x,y,z,disp(3)

  n=100
  size=10
  dx=size/n

  lambda=1.2
  mu=2.3
  slip=1
  opening=0
  width=2
  length=2
  strike=0
  dipd=90
  rake=0
  origin(1)=size/2
  origin(2)=size/2
  origin(3)=6

  do i=0,n
     do j=0,n
        x=i*dx
        y=j*dx
        z=size/2
        call okada(lambda,mu,slip,opening,origin(1),origin(2),origin(3),&
             width,length,strike,dipd,rake,x,y,z,disp)
     end do
  end do

end program fortran
